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I. INTRODUCTION 

Applications of the group theory constructions remains not broadly used yet in the kinetic theory, as it deserves 
to be [l "l In the paper, the new opportunities and advantages that the parameterization of two-particle collisions by a 
matrix from the group of rotations can provide are shown. Usually, to construct the collision integral in the 

Boltzmann kinetic equation a collision of two particles is determined by setting a direction ( 71 , 77 ^ = 1) of the 
relative velocity vector of the particles after the collision 1 

, mv + u + \v-u\n , mv + u-m\v-u\n m 

V =- 1 -—; U =- 1 -—; where m = — L , V=v-«; (1) 

1 + m 1 + m m 2 

In the paper we offer to construct the collision integral using a parameterization of a scattering by a rotation matrix 
Re 0 3 which is determined by Euler's angles (p,6, y/ [4J . 

In this case, the transformation of the velocities due to a collision in contrast to case (1), becomes a linear one: 

, mv + u + R(v-u) , mv + u — mR(v — u) 

V =--- u = ---(2) 

1 + m 1 + m 


and is representing by scattering matrix S a partitioned matrix (2x2 cells). The size of each cell obviously is 

(3x3): 





m + R 

1 -R ' 


, S = 

1 + m 

1 + m 

u 

\ 2 


m( 1 — R) 

1 + mR 



v 1 + m 

1 + 777 J 


Also in contrast to (1), that linear transformation due to det 5j = 1 provides the equality of velocity volumes 

dvdu = dv'du' 


The matrixes S(R) , R e 0 3 constitute a group: SiR^ ■ S(R 2 ) — S(R t ■ R 2 ) , S ] (R) — S(R '). 


( 4 ). 
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II. THE COLLISION INTEGRAL 


To write the collision integral the integration over relative velocity directions should be 


rein 

replaced --—> 


I* dR '' * 1 •—• * • * ■ - - j 

J ——— with integration on the invariant ( dR = dR a R = dRR 0 = dR ) measure™ over the 


An J \6n 2 .. c '~--- 0 0 

group (),: dR = d\f/d(p sin QdO, 0 < l j/ ,d ,(p< 2n , J dR = 1 6n 2 . As a result, we get the collision integral 
written in an explicitly invariant and simple form: 




1 dR 


A K 


An 


where b(y,R)= wo e M*)) 1 


is a scattering indicatrix 




v • R\ 


cosine of the scattering angle, 


F^ = foV^) = ffv)V(u) is a two-particle velocity distribution function. Taking into consideration that 
8(v' - V 0 )8(«'- M 0 ) = 8(5 ^-^ 0 ) = d(^-S(R 1 )^ 0 ) we have decomposition formula for 8 -functions: 


/(8(v-v 0 ),<5(w- h 0 )) = /(£(£ -^ 0 )) = [^-&(v 0 ,^)[8(v-v')-8(v-v 0 )], (6) 

J An 

And similar for Maxwellians: 

1 (/m 0 - v 0 ), (« - Mo )) = J(*> «o > *)fu ( v - v o )- b r i v . «o» R)/m 0 - V 0 )], (7) 

J An 

where b T (v, U 0 ,/() = J b{v,Rf¥ M (li - u 0 )du . 

At the zero temperature limit we obtain the asymptotic equality similar to property (6): 

1 {fu ( v - v o )Wm (« - «o )) r -*o = 1“ b ( v o » *]fu (v ~ K h fu (v ~ v o )]+ 0(T ), (8) 

J An 


For Maxwell’s molecules (and asymptotically in the limit T —> 0 in the general case) those two decompositions 
coincide. In the case of the Maxwell molecules the indicatrix b does not depend on a velocity and consequently 
equality (8) is satisfied at all temperatures and is a corollary of the invariance of the collision operator with constant 
collision frequency with respect to Gauss transformation 1 ' 1 : 


kT , 


e 2m ' ' l(f, x ¥) = l 


( dd_\ 

,2 m. 


f,e 


ini, “ vp 


where 



\2nkT 



m(v-v') 




, du > 


( 9 ) 



III. THE COLLISION INTEGRAL FOR DISCRETE MODELS 


The form of collision operator obtained above is very useful in an application to the discrete models^ 5 " 91 . An explicit 
expression for the collision integral suitable for the discrete models can be constructed, if we replace in collision 
integral (5) the averaging over rotations from group of rotations 0 3 with an averaging over its discrete subgroup 

consisting of finite number^ elements ( j"---> — V )• 

J 16k" K 

( 10 ) 

A K A K 

K k= i A k=l 

The expressions given by formulae (6) and (8) in this case are reduced to the finite sums: 

Ajr K / ^ v 

h (<5(v - »\> )5(« - «o)) = —lH v 0 . R k l 5 ( v - v 'o y S(v - v 0 )l , (11) 

& k =1 
and 

1 d if m 0 - v o Wm O - «o )) r _rt = ^ X 6 ( V o ’ ^ I/m ( v ~ v o )“ fu ( v “ v o )]+ °( T )- (12) 

A k =i 

Formulae (11) and (12) give an explicit representation of the collision terms in the set of equations for discrete gas 
models, which can be obtained from the initial Boltzmann kinetic equations after a substitution of distribution 
functions in them in a form of an expansion on delta-functions or on Maxwell velocity clusters. 

IV. SCATTERING VARIABLES AND AN INVARIANT ENSEMBLE OF PAIRS OF 

DISCRETE VELOCITIES. 


A state of a pair of colliding particles in the kinetic theory is usually characterized by particle velocities V and 
U. The state also can be uniquely characterized by other variables that are linear combinations of the velocities V 
and U and which transform in collisions much simpler. It is convenient to choose variables p and w as such 
scattering variables. The variable w is a center of masses velocity, and p is a momentum of the first particle in a 
center of masses frame of reference. Formulae of transformation from a variables V , U to variables p ,W and 
back look as follows: 



The matrix of scattering in the variables p, W becomes cell-diagonal. In fact, a velocity of a center of masses after 

collision does not change ( w' = W ) and the momentum p is rotating by a rotation matrix R , p = Rp. 
Consequently, the matrix of scattering in these variables has a form: 


S A = ASA.-' 



(14) 



If we will consider a model of gas in which not all possible scatterings in a pair of colliding particles 
happen, but only related to rotation matrixes Rs 0 3 (K )cO ; from some discrete subgroup 0 3 (K ) containing 

K elements, it is possible to construct an ensemble of pairs of discrete velocities, that transforms into itself by the 
collision transformations. There are several ways to construct such ensembles. We consider purely lattice way. In 
this case the velocity of a center of masses is given by an integer linear combination of the basis vectors of a lattice 

IV, ,1V 2 , IV 3 , and the momentum of the first particle is given by a linear superposition of the vectors 

m 0 IV, , m Q IV 2, m Q W 3 with integer coefficients. A set of six integer numbers l i ,l 2 ,l 3 ,k [ ,k 2 , k 3 will be 
coordinates of velocity pairs: 


p{l\,l 2 ,h)=l\m 0 w { +l 2 m 0 w 2 +l 3 m 0 w } , 
w(k u k 2 ,k 3 ) = &,iv, +k 2 w 2 +k 3 w 3 . 


vj rvj /v 2 2 

For discrete velocities in accordance to (13) one will have: 


(15) 


V a V a ih ’h’h’k\>k 2 ,k 3 ) 


^/, + k 

IV, 

+ 

^l 2 +k 2 

iv 2 + 

l 3 +k 3 




K J 


K' J 


r3 e A*; 


A i 2 ; i 3 jA; k 2 ,A)<= k81 A^2>^ 3 ? A 5 ^2’A o,±i,±2,..., 


(m ~ ) 

^ L l,R + k. 

IV, + 

^ l 2 R + k 2 

w 2 + 

^l 3 R + k 3 \ 

m a 

V “ J 


K J 


J 


( 16 ) 


l 3 R + k 3 |iv 3 , Uqp — Vp ( , l 2 , l 3 ,k^,k 2 ,k 3 ), 


The vectors IV, , IV 2 , IV 3 are basic vectors of the Bravais lattice invariant regarding a discrete subgroup of rotations 
(the point groups), 

Aiv, =^iv a T a ,(A), (17) 

k 

where T ki (r ) are unimodular matrixes with integer matrix elements. In the theory of crystals [l0 ’ n, known that there 
are 14 types of the Bravais lattices and 7 point symmetry groups of these 
lattices: S 2 ,C 2h ,D 2h ,D 3h ,D 4h , D 6h ,O h corresponding to 7 syngonies (triclinic, monoclinic, rhombic, trigonal, 
tetragonal, hexagonal, cubic). These 7 groups have 32 point subgroups. The only group not having a preferential 
direction and consequently the most appropriate one for our purposes is the O h group. That is the group of 
symmetry of a cube (isomorphic to the group of symmetry of an oclahedron) is consisting of 48 elements. Its 
subgroups without a preferential direction are: the group of proper rotations of a cube 0(24), the group of 
symmetry of a tetrahedron T d (24) , the group of proper rotations of a tetrahedron T’ (12) and the T h (24) group of 
proper rotations of a tetrahedron with an inversion added. The number of elements in these groups is indicated in 
brackets. For plane Bravais lattices there are only 10 point groups of symmetry: C l ,C 2 ,C 3 ,C 4 ,C 6 that are the 
groups of rotations around axes of the 1-st, 2-nd, 3-rd, 4-th and 6-th order (angles of rotations in the groups are 
71 Jl 

divisible by either —, or —) and the groups D ] ,D 2 ,D 3 ,D 4 ,D (i that are the groups of rotations C n with a 

reflection in an axis lying in a plane of the lattice. For a one-dimensional lattice, obviously, there are only two 
groups. Those are the identical transformation E and the group with a reflection {EE }. The fact of necessity of 
a rational masses ratio has been discovered and first lattices of discrete velocities for a mixture of gases have been 
constructed in paper [8] , paper [9] is devoted to spurious invariants in discrete models. 



Presenting velocity distribution functions f (v) by the expansions on Dirac delta-functions concentrated 
on the velocities of the discrete sets L a , 

faiy)= X/( V o)5(v-V 0 ), (18) 

v o^ a 

and two-particle distribution functions of a pair of colliding particles of the sorts (X and by the expansions on 
Dirac delta-functions concentrated on bivectors from the invariant ensembles L a p , 

fafi(£)= X/«( V »ft(“o) 5 ( V_V o) 5 (“ _ “o). (19) 


we obtain an analog of the Boltzmann equation for discrete models for gas mixtures of M components in the 
following explicit form: 


dfjy o) 

dt 


+ V 0 


dr 


= X X V X b «P ( V ’ ^ If a ( V 0 Yfn ( M 0 ) “ fa ( V 0 )fp («0 )] 

P {«ol(v 0 ,«ohi^} ^ k=] 


( 20 ) 


v 0 eL a ; = 

Where L a are sets o f discrete velocities V 0u for each mixture component depend only on the mass m a ; L a 3 are sets 
of colliding pairs (v 0a , U 0 p f G L a p ; M is number of mixture components. The definition 
{«ol(v 0 ,«o) e L a p } = Lp (v 0 ) means a set of all U 0 for which the pair | V 0 , U 0 j with the given V 0 is contained 
in the invariant ensemble L a p . 


V. EXAMPLE OF DISCRETE VELOCITY MODEL CONSTRUCTED ON THE GROUP 

OF SYMMETRY OF A CUBE O h 


Rotation on angle B around axis B is given by formula: 


„ 1-cosi? *2 sin B n , 

exp B =--— B~ +- B + 1 

B 2 B 


where the operator (matrix) B is defined by the following 


expression: 


„ df „ df 

B=Bx , Bv=Bxv- 


( 21 ) 


( 22 ) 


We will describe nontrivial rotations from the group O by vectors B ] B 13 . Rotation on the zero angle (trivial 
rotation) is the unity transformation E : 


B 0 =0, 


e x P (b q )=E 


(23) 



Symmetry of a cube includes three forth odder axises coming though centers of opposite faces of a cube (9 nontrivial 
rotations): 


{B\ 9 =±^{ 1,0,0), ±|(0,1,0), ±|(0,0,U ^(1,0,0), 7r(0,l,0), /r(0,0,l), (24) 

four third odder axises coming though opposite corners of the cube (8 nontrivial rotations): 


Ml0....,7=±^ 


+ J-,+ — 


(25) 


and six second order axises coming through centers of opposite edges of the cub (6 nontrivial rotations): 




n 


— — 

' 4l' J2 


n 


±1 


\ ( 


0 , —p= , n 


+ 11 ^ 


4i '4ij 4i 


(26) 


One can obtain the group of all transformations of symmetry of a cube O h by adding the inversion transformation 
(— E ) to the group O . All 48 elements of the group can be written making use of definitions (21)— (26) as 


{rI = { ± exp{h n ),n = 0 , 1 ,... , 23 }. 


(27) 


There are three Bravais lattices (cubic system), which are invariant under transformations from group O h . The first 
one is the simple cubic lattice F c . Sites of that lattice are situated in corners of identical cubes: 

w(k ] , k 2 , k 3 ) = k ] m>, + k 2 w 2 + & 3 , k ]t k 2 k 3 — 0, ± 1 , ± 2;:: , 

where vectors of elementary periods are 

w t =(l,0,0), w 2 = (0,1,0), w 3 = (0,0,1) . (28) 


The body centered cubic lattice T^is a lattice with sites, which are situated in vertices of cubs, and in its centers. 
Vectors of elementary periods of the lattice are three vectors from a center of the cube to any three vertices. One of 
the possible variants is: 



(~\ 

1 

1 


ri 

n 


ft 

1 

H>i = 

1 + 

2’ 

2 J 

L w> 2 = 

U’ 

2 ’ 2 J 

. = 

h’ 

2’ 2) 


(29) 


A face-centered lattice is a lattice with sites situated in comers of cubes and in centers of faces. Vectors of 

elementary periods of the lattice are three vectors from any corner of a cube to centers of three faces. One of the 
possible variants is: 


w. 


( oX * 2 =fi.o.r 

2 21 2 2 2 


( 


i. i, <2 

2 2 


(30) 



Here we will describe characteristics of discrete velocity models constructed on simple cubic lattice (28). 
Substituting elementary periods of this lattice (28) to (16) one will have a relevant discrete velocity model. For 
discrete velocities components we will have the following expressions: 


W 0 , . 

v x - - h + k ] ; 

V v =—h +k 2 ; 

v ^^Ll i+ k 3 ; 

w, 

w, 

m. 

m o , 1 

u x = - ; 

w 0 , , 

u v l 2 + k 2 , 

m o , J 

u y =-— -l 3 +k 3 ; (31) 

m 2 

m 2 

m 2 

where k i ,k 2 ,k 3 — 0,±1,... ± k m ; 

=0,±1 

It is seen from (31) that discrete velocity sets for 


components v x ,v y ,v^are the same independent sets. To have a number of different discrete velocities less then 


m , 


m, 


number of collisions (bivectors from invariant ensemble) - and- should be rational numbers: 


m n 


m n 


m L= _Pl_. = Pi . 
m 0 q ] ’ m 0 q 2 ’ 

where p t , p 2 , q 3 , q 4 are integer numbers and at list the following inequalities should be held 


(32) 


2 L +^P\’Pi; 2k m + \>q v q 2 


(33) 


To obtain all different discrete velocities of particles with mass W = — one should take indexes l and k in the 

q 

range: 

l = -l m ,...{-l m +p); k = -k m ,...k m and / = (-/,„ +p+l),.J m ; k = {k m -q),...k m 

v x = —l + k; (34) 

P 

p 

Total number of different discrete velocities v . for particles with mass m = —will be: 

q 


N. 


N. 


f p} 

yVj 

f p \ 

v^y 


= p{ 2k m + l )+q( 2l m + l )-pq> 


= {2l m + \){2k m +l). 


q< 2 k m +\, p < 21 m +1 


(35) 

q > 2 k m + 1 or p > 21 m + 1 


For illustration we here consider the specific gas mixture consisting of (Ve(20), (V0(3O) and Ar( 40 ) gases. The 
relevant discrete model has the following parameters: 

L =k m =5; m 0 = 20 ; p x = \, q x = 1; p 2 = 3, q 2 = 2 ; p 3 =2, q 3 = l; (36) 


For one-dimensional case (v = v_ =0) according to (35) the number of different velocities for each gas 
component of the mixture will be: 



(37) 


N 


^1 

lO 


= 21, N 


"2" 




= 49, N] 


V 1 7 


= 31 


The number of effective collisions (the number of bevectors from invariant ensemble with nonzero relative 
velocity, l x + 0 ) for each two gases of the mixture is 


Af t =2/.x(2t.+l) = 110 


(38) 


and 


N> 2 N 


(1' 
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"3^ 




(A 

"3 ^ 


"3^ 


'2 ^ 


"2" 
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+ N 


, N 


+ N 

- , 2 N 


, N 
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, 2 N 
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J, 
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10 

J, 


J, 


J, 


J, 


(39) 


110 >(2x21 = 42),(21+ 49 = 70),(21+ 31 = 52),(2x49 =98),(49+ 31 = 80),(2x31 = 62) 


Due to inequalities (39) this discrete model do not has spurious invariants. 

The two dimensional (v z = 0) model with l m = k m = 3 has the following characteristics: 


l m = k m = 3 ; m 0 = 20, p ] = 1, q x = 1; p 2 = 3, q 2 = 2 ■ p 3 = 2, q 2 = 1; 


(40) 


N\ 


O 

1 


= 169, N — =841, N | 


{ 2 

v 1 


= 361 ; 


N c = [(2/,„ +1) 2 ~ ijx (2k m + 1) 2 = 2352 


And calculated in the same way there are characteristics of the full three dimensional discrete model: 

L=K =2, m 0 = 20, p ] =1,^=1; p 2 =3, q 2 =2- p 3 =2,q 3 =U (41) 



rn 


f3 3 

f 2 ' 

N 


= 729, N 

- =6859 , N 



J, 

l 2 ) 

J, 


= 2197; 7V f = [(2/,„ +1 ) 3 -1 Jx (2+1 ) 3 = 15500 
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